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SECURITY  CLASSIFICATION  OF  THIS  PAOCfWh»n  Data  Bnlarad) 


CONTENTS 


I.  INTRODUCTION 

II.  THERMODYNAMIC  THEORY 

A.  The  Free  Energy 

B.  fyToelectricity 

III.  MICROSCOPIC  THEORY 

A.  Evidence  for  Disorder 

B.  Molecular  Field  Theory 

C.  Generalized  Molecular  Field  Theory 

IV.  PRESENT  STATUS  AND  RJTURE  PLANS 

A.  Experimental 

B.  Theoretical 

V.  SUMMARY 


-1- 


I.  INTRODUCTION 

The  work  described  in  this  report  was  conducted  from  1 October  1973 
to  30  September  197*+  under  Advanced  Research  Projects  Agency  Contract 
DAAB  07-7*+-C-0470,  P.M.  Raccah,  principal  investigator,  F.  Rothwarf, 
contract  monitor.  Its  purpose  was  to  improve  our  understanding  cf  the 
phenomenon  of  pyroelectricity  through  an  exploration  of  the  fundamental 
physical  processes  involved.  It  is  hoped  that  such  improved  understand- 
ing will  lead  to  practical  guidelines  for  the  selection  of  materials  for 
use  in  infrared  detectors  and  to  useful  insights  into  the  factors  limiting 
the  performance  of  devices  employing  the  pyroelectric  effect. 

The  ferroelectric  perovskites  were  chosen  as  a focus  for  this 
investigation  since  their  properties  have  been  extensively  studied,  and 
they  provide  a family  of  isomorphous  materials  in  which  the  effects  of 
varying  certain  physical  parameters  can  be  examined  systematically. 

The  pyroelectric  behavior  of  the  pervoskites  is  intrinsically  related 
to  their  ferroelectric  behavior.  Below  the  Curie  temperature  these  crystals 
are  spontaneously  polarized,  and  although  this  polarization  is  always 
screened  by  stray  charges  attracted  to  the  crystal,  any  change  in  tempera- 
ture produces  a change  in  the  polarization  of  the  crystal  which  appears  as 
a net  dipole  moment  and  which  persists  until  it  can  be  compensated  by  a 
gain  or  loss  of  stray  charge  &t  the  surfaces.  This  pyroelectric  effect  can 
be  described  by  any  theory  of  ferroelectricity  which  accounts  for  the  tem- 
perature dependence  of  the  spontaneous  polarization.  We  have  examined  tne 
simple  thermodynamic  theory  due  to  Devonshire,  and  obtained  from  it  several 
new  results  (Section  II)  which  yield  the  pyroelectric  coefficient  and  the 
figures  of  merit  for  pyroelectric  detectors  in  terms  of  the  parameters  of 
that  theory  or  in  terms  of  simple  (but  at  present  less  useful)  thermo- 
dynamic relations.  Sectior  III  describes  our  search  for  a microscopic  cal- 
culation or  the  p .rameters  appearing  in  Devonshire* s theory.  We  obtain 


the  perhaps  surprising  result  that  a suitable  generalization  of  the  mean 
field  theory  so  familiar  in  ferromagnetism  can  be  used  here  with  consi- 
derable success.  The  generalization  required  allows  the  theory  to  des- 
cribe first-order  as  well  as  second-order  transitions,  and  thus  represents 
an  interesting  contribution  to  the  general  theory  of  phase  transitions  as 
well  as  a partial  solution  to  the  present  problem.  In  justifying  the  use 
of  a mean  field  theory,  we  must  overcome  the  objection  that  the  symmetry 
of  the  high-temperature  phase  does  not  permit  the  existence  of  the  per- 
manent dipole  moments  assumed  in  the  theory.  A variety  of  experimental 
evidence  for  the  existence  of  a partially  disordered  array  of  dipoles 
in  the  pervoskites  is  presented,  including  our  own  measurements  of 
Raman  scattering  in  KNbO,.  Finally,  Section  IV  describes  additional 

■J 

experiments  now  in  progress  which  may  yield  further  evidence  for  the 
order-disoider  nature  of  the  phase  transitions  and  outlines  some  new 
theoretical  ideas  which,  although  still  speculative,  show  premise  of 
greatly  enhancing  the  power  of  the  mean  field  theory. 


| 


I 
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II.  THERMODYNAMIC  THEORY 

1 2 

Devonshire's  thermodynamic  theory  of  ferroelf  tricity  ’ can  be  used 

as  the  basis  for  a description  of  the  pyroelectric  effect  in  the  perovskites. 

Pert  A of  this  section  describes  the  free  energy  function  introduced  by 

Devonshire  and  summarizes  the  well-known  results  obtained  from  that 

function.  Part  B employs  Devonshire's  formalism  to  obtain  expressions 

for  the  pyroelectric  coefficient  and  the  figures  of  merit  relevant  to 

the  operation  of  pyroelectric  detectors.  These  latter  results,  while 

simple,  do  not  appear  to  have  been  discussed  previously. 

A.  The  Free  Energy.  If  we  restrict  our  attention  to  the  transition 

from  a ncn-polar  (paraelectric)  phase  having  inversion  symmetry  to  a polar 

(ferroelectric)  phase  and  consider  only  polarizations  P along  the  polar 

axis,  then  the  free  energy  at  zero  stress  can  be  expanded"*"  in  a power 
? 

series  in  P : 

G = £ A P2  + £ B P4  + g C P6  + . . . , (1) 

where  the  coefficients  A,  B,  C,...  are  functions  of  the  temperature. 

(Strictly  speaking,  the  free  energy  G at  non-zero  stress  is  the  elastic 
Gibbs  function.1)  The  applied  electric  field  is  given  by 

E = dG/dP  = AP  + BP3  + CP5,  + ...,  (2) 

and  the  dielectric  stiffhess  (reciprocal  of  the  susceptibility  x)  is 

1/X  = 32G/dp2  » A + 3BP2  + 5 CP4  + ...  . (3) 

For  the  susceptibility  of  the  paraelectric  (P=0)  phase  to  exhibit  a 
Curie-Weiss  behavior,  it  is  sufficient  to  assume 


-4- 


A = A (T-T  ),  (4) 

o o’ 

where  T is  the  absolute  temperature  and  T is  the  Curie -Weiss  temperature 

(For  x»l,  the  usual  Curie  coefficient  is  then  approximately  2v/kJ, 

2 

Most  discussions  of  the  thermodynamic  theory  further  assume  that  the  higher- 
order  Devonshire  coefficients,  B,C,...  are  temperature -independent  and  that 
the  expansion  of  G in  Eq.  (l)  can  be  cut  off  at  the  term  proportional  to 
P^.  Since  the  first  assumption  is  contrary  to  the  results  of  the  micro- 
scopic model  to  be  presented  in  Section  III,  we  will  avoid  its  use.  The 
second  assumption  so  greatly  simplifies  the  theory  that  we  will  adopt  it 
here,  with  the  warning  that  our  analysis  of  the  experimental  data  for 
BaTiO^  raises  some  doubt  as  to  its  validity.  We  should  also  note  that 
this  approximation  for  G requires  C > 0 in  order  to  prevent  the  minimum 
of  G from  occurring  at  infinite  polarization.  Equations  such  as  (l)-(3) 
which  are  correct  to  all  orders  in  P will  have  that  fact  indicated  by 
the  continuation  sign  (...). 

Conditions  for  the  existence  of  a spontaneous  polarization  in  the 
absence  of  an  applied  field  can  be  obtained  from  Eqs.  (l)-(4).  If  terms 

Q 

of  order  P and  higher  in  G are  neglected,  then  Eq.(2)  with  E = 0 has  the 

solutions  P = 0 and  P = P , where 

s 

P2  = - B/2C  (1  ±V  1-4AC/B2) . (5) 

s 

2 

In  order  to  have  Pg  real  we  must  have  4AC/B  < 1.  Hence  no  polar 

state  is  possible  at  temperatures  higher  than  T , where  T is  the  solution 
2 

of  4AC/B  = 1.  (For  A given  by  Eq.(4)  and  B,C  independent  of  T,  we  obtain 
2 2 

the  usual  result  , T = T + B /4A  C.)  For  the  polar  state  to  be  stable 
’ u d 'o 

(or  metastable)  we  must  have  G(P  ) be  a local  minimum  of  G,  and  hence 

s 
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s2G/ap^  = i/x  > o.  This  requires  the  choice  of  the  plus  (minus)  sign  in 

Kq.(5)  if  B is,  negative  i positive). 

We  must  now  consider  separately  the  cases  B > 0 and  B < 0.  If  B > 0 

2 

we  have  P„=0  at  T=T  , and  P is  negative  for  T < T < T . Hence  no 
s o’  s o u 

spontaneous  polarization  is  possible  above  T^,  and  for  T < Tq  the  polariza- 
tion increases  coE.tinuously  from  zero,  giving  rise  to  a second-order 
transition.  For  B < 0,  spontaneous  polarization  is  possible  for  Tq  < T < T , 

but  the  polar  state  is  metastable  as  long  as  G(P  ) > j(o)  = 0.  The  phase 

s 

transition  occurs  at  the  Curie  temperature  T^,  at  which  G(Pg)  = 0.  Using 
Eqs.(l)  and  (2)  with  G = E = 0,  we  obtain 

P 2 = - 3B/4C,  at  T = T (6) 

s'  c 

Since  the  polar  state  has  a finite  polarization  at  Tc,  the  transition  is 
first-order.  Combining  Eqs.  (5)  and  (6)  yields 

A = 3B2/l6C  , at  T = T . (7) 

(For  A given  by  Eq.(4)  and  B,C  independent  of  T,  Eq.  (7)  gives  the  usual 

result2,  T = T + 3B2/l6  AC.) 

’ c o ' o 

Equations  (6)  and  (7)  can  be  used  to  compare  the  susceptibilities 
and  Xj.  of  the  paraelectric  and  ferroelectric  phases,  respectively  . 

From  Eq.  (3)  we  have 

VXp  - A 

1/Xf  - A + 3BPS2  *50  P,. 

Direct  substitution  of  Eqs.  (6)  and  (7)  gives 

at  T = T 

c 


Xp  ■=  4 Xf  , 


(9) 
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Note  that  the  results  (5)-(9)  hold  regardless  of  the  temperatu:  dependence 

of  B and  C,  but  require  that  terms  in  G of  order  higher  than  P^  be  negligible. 

The  temperature  derivative  of  G at  constant  strain  and  polarization 
gives  the  entropy 


dG/ST  = - (|  A'F"  + £ B'P4  + r C*Pb  + ...) 


(10) 


where  the  prime  indicates  differentiation  with  respect  to  temperature. 
The  entropy  decrease  in  going  from  the  non-polar  to  the  polar  phase  is 
thus 


AS 


i A’  V c + 
* sc 


i B* 


sc 


• • • > 


(11) 


where  P is  the  value  of  P at  I = T . 
sc  sc 

Of  course,  AS  is  non-zero  only  for  a first-order  transition  (BO). 
The  latent  heat  associated  with  a first-order  transition  is  simply 


W = Tc  AS  . (12) 

B.  iVroelectricity.  The  pyroelectric  coefficient  is  just-  the 
temperature  derivative  of  the  spontaneous  polarization: 


cp  = - dPs/3T  . (13) 

Since  P satisfies  Eq.(2)  with  E = 0,  we  can  differentiate  this  equation 
s 

with  respect  to  T and  obtain 


A'P  + B'P3  + C'P5  + ... 
_ _s s s 

A + 3Bpf  + 5C  P„  . . . 

3 3 


X (A'Ps  -t  B' + C'P^  + ...),( 14) 


where  the  prime  again  indicates  differentiation  with  respect  to  temperature. 
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Variou.s  figures  of  merit  have  been  defined  for  the  operation  of 
4-6  5 

pyroelectric  detectors  . Byer  and  Roundy  have  shown  that  the  response 
of  a pyroelectric  detector  in  the  high-frequency  limit  is  proportional 
to  cp/e,  where  e = 1+4tt  x is  the  dielectric  constant.  So  for  X » 1 we 
may  define  the  response  figure  of  merit 

Mx  = <P/X  , (15) 

and  the  thermodynamic  theory  then  gives 

M,  = A'P  + B'P3  + C'P^  + ...  (16) 

1 s s s 

5 6 

It  has  also  been  shown  ’ that  the  signal  to  noise  ratio  of  a pyro- 
electric detector  is  proportional  to  co//e,  implying  a noise  figure  of 
merit 

M2  = Mx  X^  . (17) 

It  is  clear  from  Eqs.  (l4)  and  (l6)  that  any  microscopic  model  for 
ferroelectricity,  in  order  to  be  useful  in  predicting  pyroelectric 
coefficients  and  pyroelectric  detector  performance,  must  successfully 
predict  the  temperature  dependence  of  the  Devonshire  coefficients.  As  will 
be  shown  below,  the  lowest-ordeo  approximation  to  Eq.(l4),  <p  = A xP  , 
does  not  provide  an  adequate  description  for  BaTiO^  and  is  likely  to  fail 
for  other  perovskites  as  well. 

It  should  be  noted  here  that  Eq.(l4)  is  just  a special  case  of  the 
general  result 

<p/X  = (<ffl/9T)p  . u8) 

Equation  (l8)  follows  from  the  thermodynamic  relation 


1 
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,SEx  ,dl\ 

(ap l {^l 


0 , 


which  is  valid  at  constant  stress.  A general  examination  of  the  thermo- 
dynamic significance  of  the  pyroelectric  coefficient  and  figure  of  merit 
is  currently  in  progress. 


k. 


III.  MICROSCOPIC  THEORY 


Two  distinct  xypes  of  theory  have  been  developed  to  describe  the 
ferroelectric  phase  transition.  The  more  recent  a ad  generally  fruitful 
of  these  is  based  on  the  concept  of  the  lattice  dynamical  instability,  or 

y 3 

"soft  mode",  introduced  bv  Cochran'  and  Anderson  . This  theory  assumes 
the  existence  of  an  infrared-active  phonon  mode  whose  frequency  of  vi- 
bration is  small  due  to  a compensation  of  she  short-range  restoring 
forces  by  long-range  Coulomb  interactions.  By  making  reasonable  assumptions 
concerning  the  temperature  dependence  of  the  force  constants  or  by  in- 
troducing the  anhai^unicity  of  the  lattice  explicitly,  one  can  obtain  a 
soft-mode  frequency  which  vanishes  at  a characteristic  temperature  T . 

The  crystal  then  acquires  a distortion  and  polarization  of  the  same 

symmetry  as  the  soft  mode,  and  the  Lydonne -S ache - Teller  relation  predicts 

7 

a Curie-Weiss  law  for  the  dielectric  constant  . The  soft-mode  theory  thus 
establishes  a connection  between  th”  macroscopic  Devonshire  coefficients 
and  microscopic  parameters  such  as  inter-atomic  force  constants  and 
effective  charges.  Its  major  drawback  in  the  present  context  lies  in 
the  fact  that  there  do  not  exist  adequate  a priori  theories  for  th°  micro- 
scopic parameters  involved.  Instead,  one  must  infer  values  of  these  para- 
meters from  the  observed  values  of  the  macroscopic  quantities.  Hence  in 
its  present  state  of  development  the  soft  mode  theory  lacks  the  predictive 
power  needed  to  relate  the  pyroelectric  behavior  of  a crystal  to  its 
microscopic  properties. 

In  searching  for  a formalism  better  suited  to  the  problem  at  hand  we 

Q 

have  re-exa^'ned  older  theories  based  on  the  mean-field  approximation  . 

Such  theorie  de  ibe  the  ferroelectric  phase  transition  as  an  order- 
disorder  transformation  resultir  - from  the  ordering  of  an  array  of  T ermanent 
dipoles,  as  in  the  Weiss  molecular- field  theory  of  ferromagnet-sm.  This 
approach  would  not  appear  to  be  applicable  to  the  pervoskites,  where  the 
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structure  determined  for  the  paraelectric  phase  (space  group  O^1)  has 

each  atom  located  at  an  inversion  center,  thus  precluding  the  existence 

of  permanent  dipoles.  In  addition,  the  phase  transitions  predicted  by 

the  mean- field  theories  are  of  second-order,  while  many  of  the  perovsVites 

2 

exhibit  first-order  transitions  . The  purpose  of  this  section  is  to  show 
that  neither  of  these  objections  should  be  regarded  as  ruling  out  the 
possibility  of  an  adequate  description  of  ferroelectric ity  in  the  perov- 
skites  based  on  a mean-field  theory.  Part  A presents  experimental 
evidence  for  the  existence  of  disorder  in  the  paraelectric  phase  of  the 
pervoskites , indicating  that  the  atomic  sites  may  be  inversion  centers 
only  in  some  average  sense  and  that  a consistent  interpretation  can  be 
given  to  the  dipole  momeit  whii'h  appears  in  the  mean-field  theories. 

With  the  first  objection  to  this  approach  removed,  part  B outlines  the 
results  of  the  conventional  molecular-field  theory,  using  the  Lorentz 
form  for  the  local  field.  Part  C then  introduces  a generalized  form  for 
the  local  field  which  permits  the  existence  of  first-order  as  well  as 
second-order  transitions,  thus  eliminating  the  second  objection  mentioned 
above.  Although  introduced  in  an  ad  hoc  way,  this  local  field  can  be  shown 
to  represent  a reasonable  generalization  of  the  Lorentz  field  and  to 
correspond  in  a formal  way  to  a modification  of  the  molecular  field  intro- 
duced by  Onsager  in  order  to  correct  well-known  failures  of  the  Lorentz 
local  field  in  dealing  with  liquids  and  solids  containing  permanent  dipole 
moments . 

A,  Evidence  for  Disorder.  The  first  indication  of  disorder  in  the 
pervoskites  was  obtained  from  the  observation  of  diffuse  scattering  of 
electrons10  and  x-rays11  in  tetragonal  BaTiO^.  Diffuse  scattering  was  seen 

in  all  {lOOj  reciprocal  planes  except  those  passing  through  the  origin. 

12 

Comic  et  al  extended  these  observatior  to  all  four  phases  of  KNbO^  and 
found  that  as  the  temperature  was  lowered,  one  set  of  diffuse  scattering 
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planes  disappeared  at  each  phase  transition,  until  the  lowest-temperature 
(rhombohedral)  phase  exhibited  no  diffuse  scattering.  In  the  ferroelectric 
phases,  scattering  was  seen  only  in  reciprocal  planes  perpendicular  to  those 
cubic  axes  which  were  themselves  perpendicular  to  the  spontaneous  pola- 

^ i p i o 

rization.  This  behavior  was  interpreted  by  Comes  et  al  ’ , following 

14 

a model  proposed  by  Takahasi  , as  arising  from  correlated  displacements 

of  the  metal  atoms  from  the  centers  of  the  oxygen  octahedra.  The  observed 

scattering  in  the  cubic  phase  was  accounted  for  by  assuming  that  each 

metal  atom  is  displaced  along  one  of  th ; t . t equivalent  < 111  > directions 

in  which  the  displacements  along  the  chain  direction  arc  correlated  for 

distances  of  from  10  to  30  unit  cells.  Averaged  over  distances  large 

compared  to  this  correlation  length,  howev  r,  the  structure  has  cubic 

symmetry.  At  the  cubic -to -tetragonal  transition,  long-range  order  appears 

along  one  of  the  cubic  axes,  which  then  becomes  the  tetragonal  c-axis  (and 

the  direction  of  spontaneous  polarization),  while  only  the  suort-range 

order  persists  in  the  other  < 100>  directions.  The  lower -temperature 

transitions  can  then  be  described  as  the  appearance  of  long-range 

1^ 

order  along  the  remaining  cubic  axes.  It  is  readily  shown  ~ that  this 

model  accounts  for  the  observed  form  of  the  diffuse  scattering  and  for  its 

progressive  disappearance  at  each  of  the  transition  temperatures. 

More  recently,  similar  diffuse  scattering  patterns  have  been  seen  in 

15 

inelastic  neutron  scattering  experiments  on  BaTiO^,  KNbO^,  KTaO^, 

KTaxNbi  x°3  SrTiO^.  Both  neutron  scattering  and  AMR  measurements 

17 

have  shown,  however,  that  the  displacement  correlations  are  dynamic  ’ , 

% 1°  13 

rather  than  static  as  originally  proposed  by  Comes  et  al  . The  dynamic 

character  of  the  correlations  is  not  an  obstacle  to  the  use  of  an  order- 

id 

disorder  model,  however.  As  noted  by  Yamada  et  al  , the  neutron  scattering 
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results,  which  are  generally  interpreted  in  terms  of  an  overdamped  soft 
phonon,  can  be  equally  well  described  as  a tunneling  mode  i.:  which  the  atoms 
tunnel  through  a potential  barrier  between  eqi  ilibrium  positions  symme- 
trically arranged  about  the  center  of  symmetry.  In  the  latter 
picture,  the  unit-cell  dipole  moment  which  appears  in  the  molecular 
field  calculations  can  be  interpreted  as  an  instantaneous 
value  (corresponding  to  one  of  the  equilibrium  positions)  and  the  free 
energy  obtained  by  time  averaging  has  the  same  form  as  it  would  if  the 
dipole 3 were  static. 

Further  evidence  for  disorder  in  the  perovskites  has  been  obtained 

1°-21 

from  Ramar  scattering  measurements  * . No  first-order  Raman  scattering 

can  be  produced  by  a crystal  having  the  cubic  perovskite  structure,  since 

19 

each  atom  is  located  at  an  inversion  center.  Yet  Fontana  and  Lambert 
observed  broad  scattering  peaks  in  tetragonal  RaTiO^  which  snow  no  sharp 
drop  in  intensity  at  the  cubic-to-tetragonal  transition  temperature;  the 
integrated  intensity  of  these  peaks  exhibits  a temperature  dependence 

different  from  that  expected  for  either  first-  or  second-order  scattering. 

22 

Although  Barbosa  et  al  have  suggested  that  this  peculiar  temperature 

dependence  is  not  intrinsic  but  rather  due  to  an  increase  in  optical 

absorption  at  high  temperatures,  careful  measurements  of  both  absorption 

and  Raman  scattering  on  the  same,  high-quality  sample  confirm  that  the 

23 

behavior  is  in  fact  intrinsic  . 

19 

Fontana  and  Lambert  suggest  that  these  observations  could  be  accounted 
for  if  the  anomalous  features  in  the  Raman  spectra  are  interpreted  as 
disorder-indr-’ed  first-ordsr  scattering,  i.c.  one-phonon  scattering  by  modes 
which  would  not  be  Raman-active  in  a perfect  cr  al  but  which  become  so 
when  both  the  translational  and  point-group  symmetries  are  broken  by  the 
disorder.  Since  all  phonons  (not  only  those  at  the  zone  center)  can  become 
Roman  active  in  this  way,  one  anticipates  a scattering  spectrum  proportional 
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(in  first  approximation)  to  the  one-phonon  density  of  states.  We  have 
performed  Raman  scattering  measurements  on  good  quality,  single 
domain  orthorhombic  KNbO^  crystals  (provided  by  E.  Wiesendanger,  ETH 
Zurich),  and  some  of  the  results  are  shown  in  Figs.  1-3.  Figure  1 com- 
pares the  results  obtained  in  two  different  scattering  geometries, 

X(YY)X  and  X(ZZ)X.  (The  notation  ^(e^eg)kE  indicates  the  directions 
of  the  incident  and  scattered  wavevectors  k^  and  kg,  respectively,  and 
of  the  polarization  vectors  of  the  incident  and  scattered  waves  and 

—4 

e , respectively.  The  axes  are  those  of  the  high-temperature  cubic 

phase,  with  the  convention  that  the  spontaneous  polarization  P is 

s 

along  (llO).)  The  two  spectra  are  strikingly  different,  despite  the 

*4  —4 

fact  that  both  involve  the  same  momentum  transfer  k -k.  and  the  same 

s 1 

2 oa  

symmetry  component  (A^)  of  the  Raman  tensor  . In  the  X(YY)X  geometry, 
both  incident  and  scattered  photons  are  polarized  perpendicular  to  Z, 
the  direction  in  which  the  disordered  chains  are  formed,  and  the  result 
is  a typical  first-order  spectrum,  showing  sharp  peaks  against  an 
essentially  negligible  background.  The  X(ZZ)X  geometry,  on  the  other 
hand,  has  both  incident  and  scattered  photons  polarized  parallel  to  the 
z-axis,  and  the  resulting  spectrum  contains,  in  addition  to  the  peaks 
seen  in  the  X(YY)X  spectrum,  a broad,  relatively  featureless  background 
whose  intensity  is  comparable  to  that  of  the  one-phonon  peaks.  This 

background  appears  to  have  a maximum  near  150  cm  as  noted  by  Quittet 

2C  15 

ot  al  , neutron  scattering  experiments  have  located  a rather  flat 

phonon  branch,  corresponding  to  TO  modes  propagating  perpendicular  to 

the  z-axis,  near  this  energy  which  could  be  respons  ''e  for  a maximum  in 

the  density  of  states. 

Figure  2 illustrates  a second  interesting  aspect  of  our  experimental 
findings.  The  geometry  Y(XZ)Y  yields  scattering  from  modes  of  symmetry^6, 
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which  is  that  of  the  soft  mode  responsible  for  the  orthorhombic -to- 

15  

rhombohedral  transition  / In  the  Y(XZ)Y  spectrum  (Fig.  2a),  the  peak  at 
-1  20 

39  cm  has  been  attributed  to  a density  of  states  maximum  due  to  the 
15  -1 

flat  TA  branch  near  50  cm  . We  have  found,  however,  that  if  the  crystal. 

is  rotated  about  its  x-axis  (see  Fig.  2b)  this  maximum  shifts  rapidly 

toward  higher  energies,  A typical  spectrum  for  the  rotated  crystal,  labeled 

Y(XZ)Y  + R,  is  shown  in  Fig.  2a.  In  this  measurement  the  direction  of 

propagation  of  the  light  inside  the  crystal  has  been  rotated  away  from 

the  y-axis  by  less  than  12°,  but  the  position  of  the  peak  is  shifted  by 

nearly  25  cm  \ or  more  than  50$.  This  behavior  is  inconsistent  with  the 

assumption  that  the  scattering  is  due  to  the  one-phonon  density  of  states, 

and  we  believe  instead  that  it  is  produced  by  a phonon  of  nua-  -zero  wave- 

vector  whose  frequency-momentum  relation  is  highly  anisotropic.  It  has 
15 

been  shown  th-it  the  lowest-lying  .*1016  is  soft  only  for  wavevectors 

*4  -4 

in  the  xy-plane,  so  that  as  we  rotate  the  momentum  transfer  ’^-k^  (and, 
assuming  that  crystal  momentum  is  conserved,  the  wavevector  of  the  emitted 
phonon)  from  the  y-axis  toward  the  z-axis,  we  should  expect  a significant 
increase  in  the  phonon  frequency. 

23a 

It  can  be  shown,  using  arguments  similar  to  those  of  Loudon  J for 
the  case  of  a uniaxial  crystal,  that  if  the  anisotropy  of  the  short-range 
force  constants  is  small  compared  to  the  electrostatic  interactions,  the 
angular  dependence  of  the  phonon  frequency  ui  should  be  given  by 


2 

m = 


2 2 
sin  0 + 


2 2 
ip,  cos  0 


(19) 


where  0 is  the  angle  between  the  phonon  wavevector  and  the  z-axis.  We 
measured  m for  external  angles  of  incidence  (0’  in  Fig.  2b)  up  to  55°, 
corresponding  to  values  of  0 v assuming  a refractive  index  n = 2.3)  ranging 
fi'.>m  69°  to  90°.  We  also  found  u>  for  0 = 45°  from  spectra  taken  in  the 
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Y(ZX)Z  geometry.  The  results,  indicated,  by  the  triangles  and  square  in 

Fig. 3,  are  in  good  agreement  with  the  prediction  of  Eq.  (19),  shown  by 

the  solid  curve  in  Fig.  3.  A least-squares  fit  gives  = 250  cm  ^ and 

<4,  = 39  cm  \ Unfortunately  it  is  not  possible  to  observe  a mode  of 

symmetry  when  the  momentum  transfer  is  along  the  z-axis  (i.e.  9=o).  We 

have  observed,  however,  that  spectra  corresponding  to  A-^  symmetry  with 

momentum  transfer  along  z,  e.g.  Z(XX)Z,  exhibit  a broad  peak  centered  at 

250  cm  ^ and  that  this  peak  shifts  slowly  to  lower  energy  as  the  crystal 

is  rotated  about  the  x-axis.  The  angular  dependence  of  this  peak  in  the 

range  0 < 9 < 21°  was  measured  (circles  in  Fig.  3)  and  also  found  to 

fit  Eq.  (19).  This  agreement  suggests  that  we  are  in  fact  seeing  the 

same  mode  in  each  set  of  measurements,  despite  the  group  theoretical 

restriction  that  only  modes  should  be  seen  in  the  Y(XZ)Y  geometry 

and  only  A^  modes  for  Z(XX)Z.  Such  an  explanation  is  admissible  if  we 

note  that  the  group- theoretical  selection  rules  hold  only  for  phonons 

exactly  at  the  center  of  the  Brillouin  zone.  The  phonons  observed  in 

Raman  scattering  have  finite  (although  small)  v/avevectors  equal  to  the 

— ♦ — * 

mementum  difference  k -k..  At  finite  wavevector  an  infra-red  active 

s 1 
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phonon  car  have  mixed  symmetry  J , with  atomic  displacements  e ual  to  a 
linear  combination  of  displacements  transforming  according  to  each  of  the 
zone-center  symmetries  which  is  infrared  active  (in  this  case  A^,  B^,  and 
Bg).  The  results  shown  in  Fig.  3 suggest  that  this  m being  may  be  quite 
strone.  The  strongest  objection  to  our  identification  of  the  39  cm  1 line 
•as  t..e  soft  mode  is  that  the  neutron  scattering  experiments  of  Currat  et 
al"^  place  the  soft  mode  frequency  of  about  25  cm”"*"  at  room  temperature. 
This  discrepancy  could  be  caused  by  a difference  between  our  sample  and 
theirs  in  the  value  of  the  orthorhomb ic-to-rhombohedral  transition  tempe- 
rature. Unt'01  tunately,  the  transition  temperature  was  not  measured  for 


either  sample,  because  of  the  danger  of  damage  to  the  crystal  in  passing 

through  a first-order  phase  transition. 

B.  Molecular  Field  Theory.  The  evidence  presented  in  the  previous 

section  strongly  suggests  the  presence  of  disordered  dipoles  in  all  but 

24 

the  loves t -temperature  phases  of  the  pervoskites.  Lambert  and  Comes 

have  examined  the  consequences  of  assuming  that  the  phase  transitions 

in  these  materials  are  of  the  order-disorder  type  and  concluded  that  such 

a description  is  consistent  with  experiment.  In  the  following  we  will 

. 24 

re-examine  the  molecular  field  theory  (MFT)  used  by  Lambert  and  Comes 
in  order  to  show  that  it  is  in  fact  not  adequate  to  describe  the  behavior 
of  the  pervoskites  and  in  order  to  lay  the  foundations  for  the  generalized 
molecular  field  theory  to  be  presented  in  the  following  section. 

Since  the  MFT,  as  applied  to  ferromagnetism,  is  treated  in  most  ele- 
mentary texts  on  statistical  mechanics0,  we  will  simply  summarize  the 
results  obtained  in  the  case  of  ferroeicctricity. 

The  instantaneous  energy  of  orientation  of  an  electric  dipole  moment 
►i  is  given  by 

H = - p . (20) 

where  is  the  local  field  acting  on  the  dipole.  The  usual  MFT  is  ob- 
tained by  assuming  that  the  local  field  has  the  Lorentz  form 

= E + ^ P , (21) 

where  E is  the  applied  (macro icopic)  electric  field,  and  P is  the  macros- 
copic polarization.  If  we  consider  a cubic  (or  isotropic)  medium  (so  that 
P is  parallel  to  E)  and  assume  that  the  dipole  moment  can  rotate  freely, 
the  usual  statistical  a guments  for  an  array  of  N independent  dipoles  give 

P(T)  = Up  L(x) 


(22) 
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where 


L(x)  = cothx-l/x  and  x = pE^/kT  (2 3) 

(k  is  Boltzman's  constant). 

For  small  values  of  x,  L(x)  can  be  expanded  in  a power  series 

L(x)  = jx-^x3+^x5  + ...  , (2h) 

and  this  series  can  be  inverted  to  give 

x = 3y  + | y3  + y5  + ...  , (25) 

where  y = L(x)  = P/N|i.  Equation  (25)  which  relates  and  P,  can  be 
solved  for  E to  give 

E = A (T-T  ) P + BP3  + CP5  + ...  , (26) 

where 

Ao  = 3k/Np2 

B = 3A  T/N^2  (27) 

C = 99AqT/175  »V\  etc. 
and 

To  = lm/3A_  . (28) 

Since  the  electric  field  (26)  has  the  same  form  as  Eq.  (2),  we  obtain 
a >66  energy  of  the  Devonshire  form 


G = 2 Aq(T. _q)  P2  + \ BP4  + £ CP6  + . . . . (29) 

An  important  property  of  the  free  energy  (29)  is  that  the  coefficient  of 
4 ,22 

P , B = 3AQT/h  u , is  positive  for  all  finite  temperatures.  Hence  the  MFT 
predicts  a second-order  phase  transition.  This  conclusion  is  not  altered 
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if  one  replaces  the  assumption  that  the  dipole  moment  can  rotate  freely 

by  the  requirement  that  it  be  either  parallel  or  antiparallel  to  the 

2k 

applied  field  or  that  it  have  one  of  a set  of  possible  orientations 

consistent  with  the  cubic  symmetry^.  The  addition  to  Eq.  (22)  of  an 

26 

induced  polarization  dE^,  proportional  to  the  molecular  polarizability 
a,  also  has  no  effect  on  the  order  of  the  phase  transition.  The  obser- 
vation of  first-order  transitions  in  the  pervoskites  casts  serious  doubt 
on  the  validity  of  a MFT  based  on  Eq.  (21).  Since  the  theory  predicts 

a second-order  transition  which  has  no  latent  heat,  to  attempt,  as  do 
*.  24 

Lambert  and  Comes  , to  identify  the  saturation  polarization  Paflt  = Np, 

with  the  spontaneous  polarization  at  Tc  and  to  calculate  the  latent  heat 

from  Eqs.  (ll)  and  (12)  is  to  distort  the  conclusions  of  the  MFT. 

Before  discussing  our  proposed  generalization  of  the  MFr,  we  should 

note  two  other  consequences  of  the  MFT  which  are  frequently  considered 

objectionable*^.  First,  Eq.  (28)  implies  that  AqTq  = 4tt/3,  while  in 

every  ease  in  which  Aq  and  Tq  have  been  measured  in  a pervoskite,  one 

finds  a0Tq  « 4tt/3.  The  theory  thus  requires  a Lorentz  factor  in  Eq.  (2l) 

which  is  substantially  smaller  than  4tt/3,  However,  the  Lorentz  factor 

28 

varies  rapidly  with  position  in  the  perovskite  unit  cell  so  that  the 
finite  extent  of  the  dipoles  (or  the  finite  displacement  of  the  equilibrium 
positions  of  the  ions  from  the  cubic  lattice  sites)  may  produce  a subs- 
tantial redu.-.t-'  n jf  the  Lorentz  factor.  A second  objection  arises  from 
the  fact  that  if  each  unit  cell  contains  an  independent  dipole,  the 
entropy  change  produced  by  the  ordering  will  be  of  the  order  of  k per  unit 

cell,  while  experiments  indicate  a value  which  is  smaller  by  at  least  an 
27 

order  of  magnitude  . This  difficulty  can  be  met  if  the  x-ray  results  are 

inter;  -eted  to  ean  that  the  independent  r 'poles  are  actually  chains  of 

length  10-30  unit  cells,  "'e  entrop.,  change  is  then  reduced  to  approxima- 

24 

tely  the  observed  value  , although  the  conventional  MFT  continues  to 
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predict  a second-order  transition. 

C.  Generalized  Molecular  Field  Theory.  We  will  now  describe  a 
generalized  molecular  field  theory  (GMFT),  which,  by  permitting  a first- 
order  phase  transition,  overcomes  the  main  objection  raised  to  the 
conventional  MFT.  The  GMFT  is  obtained  by  simply  replacing  the  Lorentz 
local  field  (2l)  with  a local  field  of  the  form 

E^  = E + IP  + X'P3  + I'T5  + ...  . (29) 

We  have  used  the  cubic  symmetry  of  the  crystal  and  the  assumption 
that  the  polarization  occurs  along  a cubic  axis  in  order  to  suppress  the 
tensor  nature  of  X,X',X'',...  and  permit  them  to  be  written  as  simple 
coefficients.  With  this  modification  of  the  local  field,  Eq.  (25) 
now  gives 

E = Aq(T-To)P  + Bq(T-T1)P3  + C^T-T^P5  + ...  , (30) 

where 

Aq  = 3k/%2  , To  = X/Ao 

Bo  = 3Aq/5nV  , T1  = X’/B0  (31) 

Co  = 99  Aq/175  Ai4  , T2  = X"/Co 

and  so  forth.  This  leads  to  the  free  energy 

G = \ Ao(T-Tq)P2  + £ Bq(T-T1)P4  + | Cq  (T-T2)P6  + ...  . (32) 

4 

One  sees  immediately  that  for  T^  > Tq  the  coefficient  of  P , 

B = Bo  (T-Tj)  , (33) 

will  be  negative  in  the  vicinity  of  Tq,  and  a first-order  transition 
is  possible.  The  temperature  dependence  of  B has  been  measured  only  in 
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BaTiO^.  Drougard  et  al  studied  the  polarization  dependence  of  l/x  in 
the  paraelectric  phase  and  found  that  B has  the  predicted  form  (33)  with 
Bq  = 1.8  x 10  ^ cgs  and  448  K.  Also,  Huihregtse  and  Ynung^ 
determined  B from  measurements  of  the  susceptibility  and  spontaneous 
polarization  in  the  ferroelectric  phase  (just  above  the  tetragonal- 

orthorhombic  transition)  and  found  a value  in  good  agreement  with  Eq.(33) 

29 

using  the  parameters  given  by  Drougard  et  al  . More  recently,  the  work 

O] 

of  Matsuda  and  AbeJ  has  failed  to  confirm  Eq.  (33) > but  these  measure- 
ments yield  B only  indirectly,  via  a calculation  which  makes  the  dubious 
assumption  that  C is  temperature -independent. 

Before  entering  into  a detailed  examination  of  the  consequences  of 
the  free  energy  (32)  predicted  by  che  GMFT,  we  will  discuss  briefly  the 
question  of  the  local  field  in  order  to  justify  the  form  proposed  in 
Eq.  (29).  It  has  long  been  known^2  that  the  Lorentz  local  field  (21) 
ir  not  suitable  for  the  description  of  liquids  composed  of  polar  molecules 
(i.e.  molecules  having  permanent  electric  dipole  moments).  The  reason  is 
that  if  the  observed  densities  and  molecular  dipole  moments  are  substituted 
in  Eqs.  (27)  and  (28)  the  MFf  predicts  that  one  should  observe  ferroelec- 
tric ordering  in  many  common  liquids  at  room  temperature.  (Water,  for 
02 

example J , would  order  at  Tq  = 1140  K.) 

Among  the  numerous  attempts  to  avoid  this  difficulty  in  the  MFT,  one 

33 

finds  a particularly  interesting  and  elegant  approach  due  to  Onsager  . 
Consider  a medium  of  dielectric  constant  e in  an  applied  electrid  field  E. 
The  field  inside  a spherical  cavity  of  radius  a containing  a single  dipole 

For  m-  = 0 the  applied  field 
The  presence  of  the  dipole 


— 34 

of  moment  p,  can  be  calculated  as  follows 

3e  - 

produces  a field  in  the  cavity,  G = — E. 
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gives  rise  to  an  additional  field  R = ^~r  ihis  "reaction  field" 

2e+1  a3 

is  the  difference  between  the  field  which  the  dipole  would  produce  in 

vacuum  and  that  actually  produced  in  the  medium.  The  field  in  the  cavity 

is  then  the  field  of  the  dipole  itself  plus  the  sum  G+R.  If  were  an 

induced  dipole  it  wo^xld  be  related  to  the  macroscopic  polarization  by 
— 4 3 

P = |±/V,  where  V = — rra  , and  one  would  obtain  the  Loren tz  result, 

JO 

G+R  = E + 4tt/3  P.  For  a permanent  dipole,  however,  Onsager  pointed  out 

— ♦ — • 

that  since  R is  always  parallel  to  i1,,  the  only  field  tending  to  orient 

— • — ♦ ^ 

(j,  parallel  to  P is  the  cavity  field  G.  Thus  the  local  field  appearing 
in  the  orientation  energy  (20)  should  be  = G or 


E + 


4tt 

2e+l 


(34) 


where  we  have  used  P = XE  = E.  If  e » 1,  we  have  4tt/(2c+1)  ~ l/2x> 
and  the  Devonshire  result  for  the  dielectric  stiffness  (3)  gives 


E^  = E + | AP  + | BP3  + | CP5  . . . . (35) 

Thus  in  a ferroelectric  the  Onsager  local  field  is  not  linear  in 
the  polarization  but  contains  terms  of  all  orders  in  P.  While  Eq.  (35) 
can  be  regarded  as  a special  case  of  (29),  several  points  should  be  kept 
in  mind.  In  our  development  we  have  assumed  implicitly  that  the  coeffi- 
cients X,X',X",...  are  temperature-independent.  This  is  net  th°  case  for 
the  Onsager  field  (35),  and  the  consequences  are  interesting.  Although 

Onsager  believed  that  Eq.(35)  uaea  not  lead  to  a phase  transition  since 

33  35 

the  dielectric  constant  cannot  diverge  , it  was  later  shown  that  there 

is  in  fact  a first-order  transition  but  without  a Curie-Weiss  behavior  of 

35 

the  dielectric  constant.  Pirenne  has  compared  the  Lorentz  and  Onsager 
forms  of  the  local  field  in  terms  of  the  degree  of  correlation  among  the 
dipoles  which  they  imply.  The  Lorentz  form  assumes  no  correlation:  the 


1 
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dipoles  of  the  medium  do  not  respond  to  the  orientation  of  the  test  dipole, 

0n  g ^ 

so  that  the  reaction  field  is  always  equal  to  — P.  On  the  other 

hand,  the  Ousager  formula  assumes  "complete"  correlation:  the  medium  is 

always  polarized  in  a way  consistent  with  the  instantaneous  orientation 

of  the  test  dipole,  so  that  R is  always  parallel  to  p..  Since  processes 

involving  the  rotation  of  dipoles  or  the  tunneling  of  ions  are  generally 

much  slower  than  electronic  processes,  one  might  expect  the  real  situation 

q-i 

to  lie  somewhere  between  these  two  extremes.  Finally,  Frflhlich  has 
pointed  out  that  in  deriving  either  form  of  the  local  field  one  neglects 
short-range  forces  end  considers  only  long-range,  dipole-dipole  inter- 
actions.  Yet  lattice  dynamical  calculations’3  indicate  that  short-range 
forces  play  a significant  role  in  the  correlations  which  lead  to  the 
formation  of  chaint . They  are  also  likely  to  affect  the  interaction 
between  chains,  causing  a departure  from  the  simple  dipole-dipole  inter- 
action assumed  in  deriving  the  Onsager  result  (3*0.  In  summary,  the 
proposed  form  (29)  for  E^  resembles  the  Onsager  result  in  that  its  non- 
linear dependence  on  P permits  the  existence  of  a first-order  transition, 
while  the  introduction  of  the  X,X',A" , ...  ...  adjustable  parameters  permits 

the  theory  to  take  into  account  the  short-range  forces  in  an  empirical  way. 
Calculations  are  now  under  way  in  an  attempt  to  relate  the  expansion  co- 
efficients X,X',X",  ...  to  the  microscopic  quantities  which  determine  the 
dynamic  behavior  of  the  lattice. 

In  order  to  compare  the  predictions  of  the  free  energy  G (32)  with 
experiment,  we  note  first  that  the  coefficients  Aq,  Bq,  CqJ  ...  given  by 
Eq.  (31)  depend  on  just  two  parameters,  N and  n.  If  we  ignore  te-  ms  in 
G of  order  higher  than  P^,  we  have  only  three  additional  parameters,  the 

temperatures  Tq,  T^,  and  Tp.  For  BaTiO^,  four  of  these  quantities  have 

29  37 

been  determined  directly  from  experiment:  ’ 
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A = 7.6  x 10‘5  cgs  , T = 383  K 
o o 

-14  (36) 

B = 1.8  X 10  cgs  , T = 448  K 

O 1 

2 

We  obtain  from  these  the  number  of  independent  dipoles,  N = = 

20  -3  1- 

4.6  x 10  cm  J and  the  saturation  polarization  Pga+  = Np  = (3Ao/5Bq)2 

2 22 
= 17  pC/cm  . The  observed  density  of  BaTiO^  gives  Ko=  1.5  x 10  unit 

O 

cell  "-/cm0,  so  that  the  calculated  value  of  W corresponds  to  a chain 

length  L = Nq/N  = 33  unit  cells,  which  is  quite  consistent  with  the 

12  13 

observed  width  of  the  diffuse  x-ray  scattering  streaks  ’ . Experimental 

2 

estimates  of  the  saturation  polarization  lie  in  tne  range  16-26  pC/cm  , 

29 

depending  on  the  quality  of  the  sample  . Our  calculated  value  of  P ^ 
is  somewhat  low,  especially  sin-.e  the  data  (36)  are  for  a sample  wish 
a relatively  large  spontaneous  polarization,  but  the  agreement  with 
experiment  is  still  satisfactory.  In  principle,  one  more  piece  of  ex- 
perimental data  would  be  required  for  each  of  the  additional  characteristic 
temperatures  T,,,^,  ...  we  wish  to  determine.  We  can  proceed  in  an 

approximate  way,  however,  if  we  note  that  the  observed  values  of  P near 

s 

the  Curie  point  are  quite  close  to  our  calculated  P and  that,  as  noted 

Sells 

qU 

elsewhere"5  , the  approach  of  P to  its  saturation  value  as  the  temperature 

o 

is  lowered  is  much  more  rapid  than  is  predicted  by  Eq.(5)  with  constant 
values  of  B and  C.  These  observations  suggest  that  if  we  ignore  terms 
in  G of  order  higher  than  P^,  we  should  approximate  P at  T = T by  P , , 
so  that  Eq.(6)  gives 

a = - 44/35  , at  T = T,  , (37) 

where  a = (T-T^J/fT-Tg).  Evaluating  2q.(5)  at  the  Curie  point  gives 

t/o  = 21/176  , at  T - Tc  , (38) 

where  t = (T-Tq)/(T-T^)  . Equation  (38)  yields 
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T 

c 


176  T - 21  (7  T, 
0 c 1 

176  - 21  a 
c 


(39) 


where  the  subscript  c indicates  T = Tq.  Substituting  the  approximation 
(37)  and  the  experimental  data  (36)  in  Ea.  (39)  gives 


Tq  = (20T  +3T  )/23  = 391.5  K , (40) 

29  37 

which  is  in  excellent  agreement  with  the  observed  Curie  temperature  , 
T - 392  K.  We  must  point  out,  however,  that  the  approximation  (37) 
also  leads  to 


T2  = Tc  + TiV°c  = 346.5  K. 


(41) 


When  the  free  energy  (with  B < 0)  is  truncated  at  the  PD  term,  it  no 

longer  has  a stable  minimum  at  P = P if  C = C (T-T0)  vanishes,  (see 

Section  II).  Hence  the  temperature  is  a lower  bound  for  the  validity 

uf  the  assumption  that  the  higher-order  terms  are  negligible,  and  we 

find  that  results  depending  on  this  assumption,  such  as  Eqs.  (5)  and  (l4), 

can  be  correct  in  at  most  the  limited  range  TC"T  5 45  K.  In  fact,  even 

close  to  the  Curie  point  there  is  evidence  for  the  importance  of  higher- 

crder  terms:  Eq.  (9),  relating  the  susceptibilities  of  the  paraelectric 

1 39 

and  ferroelectric  phases  x^  = 4xf,  is  not  well  satisfied  . If,  despite 
these  warnings,  we  continue  to  neglect  the  higher-order  terms,  vre  obtain 
from  Eqs.  (ll)  and  (16)  the  entropy  change  at  T 


o sc 

and  the  response  figure  of  merit, 


F2  + 5 B p1*  +tC  P~  , 

* cn  4 O SC  D O SC 


(42) 


M,  = A P + BP3  + CP 
1 c s os  OS 


(43) 
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The  predictions  of  Eqs.  (42)  and  (43) , with  P = P = P , are  compared 

S SC  S£Lu 

2 27  40 

with  experiment  ’ ’ in  Table  I,  which  summarizes  the  experimental 

tests  of  the  GMFT.  The  overall  agreement  for  BaTiC^  is  quite  good, 
despite  possible  errors  arising  from  the  truncation  of  the  power  series 
for  G.  We  are  currently  attempting  to  estimate  the  magnitudes  of  such 
errors  and  to  find  methods  of  calculation  which  avoid  them. 

As  noted  in  Section  I,  the  GMFT  leads  to  contributions  to  cp  and 
from  the  temperature  dependence  of  B and  C,  and  these  are  not  negligible. 
For  example,  if  B and  C are  assumed  temperature- independent,  one  obtains 
M,  = A ? = 3.9,  which  is  less  than  half  the  result  given  by  Eq.  (43)  and 

X o s 

substantially  less  than  the  experimental  figure  of  merit  (see  Table  I). 

It  is  also  worth  pointing  out  the  important  role  of  the  local  field  in 

determining  M^.  From  Eq.(43),  together  with  Eq.(2)  and  the  observation 

that  E=0  for  P=P  , we  obtain 
s’ 

~ (XPg  + I’Pj  + + ...  ) , (44) 

where  we  can  recognize  that  the  term  in  parentheses  is  equal  to  E^  for 

E=0.  The  first  term  of  (44),  M,  = A T P /T,  was  obtained  from  the  con- 

1 o o s' 

ventional  MFT  by  van  der  Ziel^,  who  concluded  that  the  factor  Tq/T  implied 
that  increased  Curie  temperatures  would  result  in  improved  room-temperature 
detectivity.  In  fact,  the  correct  temperature  dependence  of  is  given  by 
Eq.  (43):  M,  depends  on  temperature  only  through  P . Hie  l/T  dependence 
given  by  (44),  while  correct,  is  misleading  since  it  reflects  the  fact  that 
the  orientational  contribution  to  x vanishes  as  the  polarization  saturates. 
Detector  performance  is  then  limited  by  electronic  contributions  to  X which 


are  not  included  in  the  model. 
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IV.  PRESENT  STATUS  AMD  HJTURE  FLAMS 

A,  Experimental  - In  addition  to  our  continuing  investigation  of 

Raman  scattering  in  KNbO^,  measurements  are  now  in  progress  on  BaTiO^ 

and  LiTaO^.  In  KNbO^  we  have  obtained  spectra  exhibiting  line  shapes 

characteristic  of  resonant  interference,  similar  to  the  results  of  Rousseau 
4l 

and  Porto  for  BaTiO^.  We  are  currently  in  the  process  of  fitting  and 

42 

interpreting  these  spectra  with  the  aid  of  the  theory  developed  by  Fano 
for  Auger  processes.  An  interesting  preliminary  observation  is  that 
the  interference  seems  to  involve  a sharp  line  of  (predominantly)  A^ 
symmetry  at  about  195  cm-1  and  the  same  broad,  strongly  orientation  depen- 
dent peak  discussed  in  Section  III  A.  If  our  interpretation  of  this 

latter  feature  is  correct,  we  are  seeing  interference  between  scattering 

43 

processes  involving  two  different  one-phonon  final  states  , not  a one- 

4l 

phonon  state  and  a two-phonon  state  as  proposed  by  Rousseau  and  Porto  . 

We  have  also  begun  a second  series  of  experiments  designed  to 
explore  the  nature  of  the  disordex'  in  the  perovskites.  We  plan  to  measure 
optical  second  harmonic  gene,  at  ion  (SHG)^14  as  a function  of  temperature 
near  the  cubic-to-tetragonal  transition.  Since  the  cubic  perovskite 
structure  has  inversion  symmetry,  a perfectly  ordered  crystal  having  this 

44 

structure  would  produce  no  SHG  . Thus  observation  of  SHG  in  the  high- 
temperature  phase  would  add  to  the  evidence  favoring  the  existence  of 
disorder,  and  its  spatial  distributio.  would  contain  information  concerning 
any  correlations  of  the  atomic  displacements.  SHG  measurements  below  the 
transition  temperature  will  be  of  value  also,  since  only  a limited  amount 

45 

of  data  on  the  temperature  dependence  of  the  SHG  coefficients  is  available  . 

46 

Simple  theoretical  models  suggest  that  these  coefficients  should  be 

directly  proportions]  to  the  antisymmetric  part  of  crystalline  potentional 

and  hence  to  P . This  behavior  has  been  found  in  BaTiCL  and  LiTaCL  but 
s 3 3 
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has  not  been  verified  in  general  . The  linear  electro-optic  coefficient, 
another  nonlinear  optical  coefficient  similar  to  the  SHG  coefficient 

47 

but  describing  the  mixing  of  photons  one  of  which  has  :ero  frequency  1 , 
can  also  be  shown  to  be  proportional  to  Ps  in  the  models  employed  in 
Ref.  46.  If  we  keep  only  the  first  term  in  the  expansion  (l4)  for  the 
pyroelectric  coefficient  cp,  we  obtain  the  proportionality  between  the 

48 

pyroelectric  and  linear  electro-optic  coefficients  proposed  by  Soref  . 

Equation  (l4)  shows  that  this  proportionality  should  fail  when  Pg 

approaches  its  saturation  value,  and  higher-order  terms  in  cp  become 

important.  We  intend  to  pursue,  both  theoretically  and  experimental ly, 

these  relations  between  the  pyroelectric  and  nonlinear  optical  coefficients. 

B.  Theoretical.  The  G'FP  proposed  in  Sec  cion  III  lends  itself  to 

extension  and  improvement  in  two  distinct  areas  defined  by  the  structure 

of  the  theory.  As  one  can  see  from  Eq.  (3l),  the  coefficients  Ao,Bo,Co,... 

depend  only  on  the  parameters  N and  p.  and  our  choice  of  the  function 

L(x)  = cothx  - l/x.  On  the  other  hand,  the  temperatures  ... 

depend  on  the  local  field  as  well  as  on  the  coefficients  A ,B  ,C  ,...  . 

o’  o’  o’ 

It  therefore  makes  sense  to  consider  separately  the  problem  of  determining 
the  temperatures  and  that  of  determining  the  coefficients  since,  while 
not  inde pendent , the  two  problems  are  couj led  in  a particularly  simple 
way.  Our  ultimate  goal  is  to  relate  all  these  parameters  to  microscopic 
quantities  which  can  be  obtained  from  either  experimental  or  theoretical 
descriptions  of  the  lattice  dynamics  of  the  pervoskites. 

Calculation  of  the  coefficients  in  Eq.(3l)  can  itself  be  divided  into 
two  subproblems.  The  parameters  N and  p can  be  obtained  from  lattice 
dynamics  calculations  since  the  former  is  related  to  the  correlation  length 
characteristic  of  the  chains'5  , and  the  latter  can  be  calculated  from  the 
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effective  ionic  charges  used  in  the  calculation  and  the  atomic  dis- 
placements observed  below  the  phase  transition.  Since  lattice  dynamics 
calculations  generally  require  the  input  of  a substantial  number  of 
experimentally  determined  parameters,  we  do  not  regard  this  approach 
as  inherently  superior  (or  more  fundamental)  th^i  our  pro.  dure  in 
Section  III  of  taking  uhe  coefficients  Aq  and  Bq  directly  from  ex- 
periment. The  higher-order  coefficients  Cp,  etc.  then  uepend  ' n the 
choice  of  L(x),  which  introd’ces  our  secc"d  subproblem.  The  present 
V' r\  has  used  L(x)  = cotw  - l/x,  whJch  corresponds  to  the  classical 
or  "infinite- spin"  limit  in  which  the  dipole  is  free  to  assume  any 
arbitrary  orientation  with  respect  to  the  applied  field.  The  assumption 
of  Comes  et  al  ’ J that  in  the  cubic  phase  the  unit  cell  dipole  orien- 
tations Eire  equally  divided  among  the  < 111  > directions  leads  to  a 
correspondence  with  the  spin  \ case,  in  which  L(x)  - tanhx.  We  have 
found,  however,  that  this  assumption  for  L(x)  lead3  to  poorer  agreement 
with  experiment  than  the  classical  one  and  that  the  gross  features  of 
the  diffuse  x-ray  scattering  can  be  explained  using  other  assumptions 
for  the  possible  orientations  of  the  unit  cell  dipoles.  We  are  now 
exploring  systematically  the  consequences  of  assuming  possible  dipole 
orientations  corresponding  to  higher  spin  values.  It  is  important  to  note, 
however,  that  van  der  Ziel^,  using  a conventional  MFT  and  the  assumption 
that  the  dipoles  may  be  parallel  or  anti-parallel  to  each  of  the  cubic 
axes,  has  obtained  a partition  function,  and  hence  an  expression  for  L(x), 
whJ(v»  doe3  not  correspond  to  that  obtained  for  any  value  of  the  spin. 

Thus  th  catalogue  of  spins  J = ^,1,...,  ® does  not  exhaust  the  list  of 
acceptable  partition  functions,  and  we  are  now  seeking  both  a method  of 
parametrising  * complete  set  of  such  partition  functions  and  a physical 
argument,  perhaps  .ased  on  the  d tails  of  the  diffuse  x-ray  scattering, 
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which  will  determine  unambiguously  the  correct  choice. 

Calculation  of  the  temperatures  in  Eq.(3l)  involves  the  local 

field  and  here  we  are  exploring  an  approach  which,  while  still 

speculative,  seems  at  present  to  be  very  promising.  The  local  field  (29) 

is  given  to  us  as  a power  series  in  P,  where  the  coefficients  must  be 

determined  from  experiment.  If  this  expansion  is  regarded  as  merely 

an  approximation,  it  is  reasonable  to  examine  other  fractional  forms  for 

E^  which  yield  Eq.  (29)  as  an  expansion  for  small  values  of  P.  One  set 

, 1+9 

of  possibilities  is  provided  by  the  theory  of  Pade  approximants 
Instead  of  Eq.  (29)  we  can  write 


E^  = E + [m,n]  , (44) 

where  [m,n]  = p^/q^,  the  (m,n)  Fade  approximant  for  the  power  series 
XP  + X'P3  + + ...  , is  the  ratio  of  polynomials  p^  and  which  are 

of  order  m and  n in  P,  respectively.  The  approximant  [m,n]  agrees  with 
the  series  up  to  the  term  of  order  I,r‘+n.  Ihe  lowest-order  non- trivial 
approximant  is  [2,2]  which  gives 


Et  = E + 


XP 


(45) 


,4 


1-(X'/X)P 

Ibis  agrees  with  Eq.  (29)  up  to  the  term  in  P^  (which  happens  to  vanish) 
and  predicts 

X"  = X,2/X  . (46) 


2 

Using  the  experimental  data  (36)  andEq.(3l)  we  obtain  X"  = (BqT^)  /AqTq  = 

-22 

2.23  x 10  . This  is  to  be  compared  with  the  result  of  the  GMFT 

with  L(x)  = cothx  - l/x,  which  gives  Cq  = 11  Bq2/7Ao,  T^  = 346.5  K (Eq.(4l)), 

-22 

and  X"  = = 2.32  x 10  . This  agreement  of  better  than  5$  is  strong 

evidence  that  the  coefficients  of’  P in  the  expansion  of  E^(29)  are  not 
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independent  but  are  related  as  implied  by  Eq.(45).  Equation  (45)  poses 
an  important  question,  however,  which  is  to  justify  on  physical  grounds 
the  presence  of  a pole  in  E^  and  to  determine  its  location.  The  smooth 
dependence  of  P on  temperature  would  rule  out  the  possibility  of  a 
pole  for  P < Nn  = P If  we  make  the  assumption,  as  yet  unjustified, 
that  the  pole  occurs  at  P = P we  obtain 

Paat  - W)4  W) 

2 

The  experimental  data  (36)  then  give  P ^ = 20  ^/cm  in  reasonable 
agreement  with  the  prediction  of  the  GMFT  and  with  experiment  (Table  i). 
The  success  of  Eqs.  (46)  and  (47)  suggests  that  Eq.  (45)  contains  an 
important  insight  into  the  behavior  of  the  local  field,  but  this  result 
remains  to  be  established  cn  a sound  physical  basis. 
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V.  SUMMARY 

The  work  described  in  this  report  was  undertaken  with  the  goal 
of  improving  our  basic  knowledge  of  the  phenomenon  of  pyroelectricity- 
in  a way  which  would  lead  to  practical  guidelines  for  the  selection 
of  infrared  detector  materials  and  to  an  understanding  of  the  factors 
limiting  the  performance  of  pyroelectric  detectors.  The  pervoskite 
family,  an  extensively  studied  group  of  isomorphous  materials,  was 
chosen  as  the  object  of  this  investigation  both  for  its  simple  crystal 
structure  arid  for  the  range  of  physical  properties  exhibited  by  its 
members. 

We  have  developed  a theoretical  treatment  of  pyroelectricity  in 
the  pervoskites  which  is  based  on  the  connection  between  pyroelectricity 
and  ferroelectricity  in  these  materials  (Section  II).  With  the  aid  of 
this  theory  it  was  shown  that  the  figure  of  merit  commonly  used  to 
characterize  pyroelectric  detectors  is  not  merely  an  empirical  quantity 
but  one  of  fundamental  thermodynamic  significance.  This  opens  the  way 
to  an  analysis  which  may  lead  to  a priori  limitations  on  detector 
efficiency.  The  theory  also  permits  the  figure  of  merit  to  be  calculated 
using  measurements  of  the  static  dielectric  constant  above  the  Curie 
temperature  instead  of  the  more  difficult  direct  measurements  of  the 
pyroelectric  coefficient  and  dielectric  constant  below  the  Curie  point. 

We  have  also  devised  a simpl°  microscopic  model  (Section  III), 
a generalization  of  the  molecular  field  theory  so  familiar  in  ferromagnetism, 
which  accurately  describes  the  ferroelectric  phase  transition  using  only 
a small  number  of  experimentally  accessible  parameters.  This  model  not 
only  predicts  the  ferroelectric  and  pyroelectric  behavior  of  interest  in 
the  present  work  but  is  broadly  applicable  to  structural  and  magnetic 
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phase  transitions  of  all  types.  It  is  simple  to  use  and  likely  to  be 
of  value  in  many  problems  involving  phase  transitions.  Ihe  model  has 
basic  features  which  are  quite  novel,  and  it  represents  interesting 
contribution  to  the  general  theory  of  phase  transitions. 

Finally,  we  have  performed  Raman  3cattering  measurements  on 
KNbO^  which  reveal  unexpected,  anomalous  features  that  caj.iot  be 
explained  in  term'?  of  the  conventional  view  of  the  structural  phase 
transition  but  wi  ich  are  fully  consistent  with  our  generalized  molecular 
field  theory. 
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Table  I.  Comparison  of  predictions  of  the  CM  FT  with  experimental 
results  from  Refs.  2,  27  and  40. 


Psat 
Tc  (S) 

AS  (k/unit  cell) 
W (cal. /mole) 

M1  (cgs) 


Theory 

17 

391.5 

0.03 

24 

8-5  <VPsat> 


Experiment 

16-26 

392 

0.011-0.06 

9-47 

io  (373  y ) 
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Figure  Captions 

Fig.  1.  Comparison  of  KNbO^  Reman  spectra  for  the  scattering  geometries 

X(YY)X  and  X(ZZ)X,  showing  disorder -induced  one-phonon  scattering 
for  incident  and  scattered  photons  polarized  along  the  z-axis 
(direction  of  chain  formation). 

Fig.  2.  (a)  Raman  spectra  of  KNbO^  for  the  scattering  geometry  Y(XZ)Y 

and  for  the  same  geometry  but  with  the  crystal  rotated  by  9' =27° 
about  the  x-axis  . The  latter  spectrum  is  labeled  Y(XZ)Y  + R. 

(b)  Scattering  geometries  for  the  spectra  in  (a),  showing  the 
internal  and  external  propagation  angles  9 and  9',  respectively. 

Fig.  3.  Angular  dependences  of  the  Raman  scattering  peaks  observed  at 
39  om  ^ in  the  Y(XZ)Y  geometry  and  at  250  cm  ^ in  the  Z(XX)Z 
geometry.  The  solid  line  is  the  prediction  of  Eq.  (ly). 
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